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Quark propagators of the unquenched Kogut-Susskind(KS) fermion obtained from the gauge 
configurations of the MILC collaboration are measured after Landau gauge fixing and using the 
Staple+Naik action. Presence of the qq condensates and A 2 condensates in the dynamical mass M(q) 
and the quark wave function renormalization Z^(q 2 ) are investigated. We obtain the correlation 
of the renormalization factor of the running coupling taken at /i ~ 6GeV and that of the quark 
wave function renormalization Z^(q 2 ) of the Staple+Naik action. The mass function M(q) is finite 
at q = and its chiral limit is ~ 0.38GeV. We compared the results corrected by the scale of the 
vertex renormalization and the tadpole renormalization with the corresponding values obtained by 
the Asqtad action without renormalization and observed good agreement. Implication of infrared 
finite Zi((j) = 1/Z^(q 2 ) to the Kugo-Ojima confinement criterion is discussed. 
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I. INTRODUCTION 

The mechanism of dynamical chiral symmetry break- 
ing and confinement is one of the most fundamental prob- 
lem of hadron physics. The propagator of dynamical 
quarks in the infrared region provides information on dy- 
namical chiral symmetry breaking and confinement. In 
the previous paper 1], we measured gluon propagators 
and ghost propagators of unquenched gauge configura- 
tions obtained with quark actions of Wilson fermions 
(JLQCD/CP-PACS) and those of Kogut-Susskind(KS) 
fermions (MILC) in Landau gauge and observed that the 
configurations of the KS fermion are closer to the chiral 
limit than those of Wilson fermions. 

In the analysis of running coupling obtained from the 
gluon propagator and the ghost propagator, with use of 
the operator product expansion (OPE) of the Green func- 
tion, we observed possible contribution of the quark con- 
densates, A 2 condensates and gluon condensates in the 
configurations of the KS fermionQ]. The quark propa- 
gator of quenched KS fermion was already measured in 
2| , and possible contribution of these condensates are re- 
ported. Unquenched KS fermion propagator of 20 3 x 64 
lattice (MILCc) was measured in @, but to distinguish 
the gluon condensates and the quark condensates, it is 
desirable to measure the quark propagator of larger lat- 
tice (MILC/) and to compare with data of MILC C . We 
measure quark propagator of gauge configuration pro- 
duced by using the Asqtad action 1) MILC C 20 3 x 64, 
/3 = 6.76 and 6.83 and 2) MILC/ 28 3 x 96, (3 = 7.09 and 



7.11. @ The paper is organized as follows. In sect. 2 we 
present kinematics of the staggered fermion on the lat- 
tice and in sect. 3 we show the renormalization effects and 
in sect. 4 numerical results of the mass function and the 
quark propagator are presented. Finally in sect. 5 we give 
conclusion and discussion. In the appendix we show the 
algorithm of the inversion to get the quark propagator. 



II. KINEMATICS OF THE KS FERMION 

The methods of writing the KS fermionic action on 
the lattice and deriving its propagator are presented in 
0. The KS fermion contains 16 flavor degrees of freedom 
and it is necessary to reduce the extra degrees of freedom 
by taking the 4th root of the fermion determinant. The 
MILC collaboration improved the fermionic action by fat- 
tening the gauge links, including Naik term and multi- 
plying the tadpole renormalization factor so that the con- 
tinuum limit can be approached on relatively small size 
of the lattice^, 0. Although there are no rigorous jus- 
tification, there are indications that taking the 4th root 
does not yield serious problems. 

In 3] , a formulation of the KS fermion propagator cal- 
culation was presented and the numerical calculation was 
performed on MILC C (20 3 x 64 lattices) data0. In these 
papers, the momentum of quarks on the hypercubic lat- 
tice is defined by 



L 



where n„ = 1, 2, • • • L/4 and the 16 flavor degrees of free- 
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In the Asqtad action, the link variable is modified by 
fattening 

V 

+ E^ 5 ^) + E (*))] (3) 



(3) 

where Sp,J is the staple contribution 

S$(x) = U v {x)U^{x + v)Ut{x + A) + (4) 

S^ 5 * 1 and S^ 7 - 1 are 5-link and 7-link contribution, respec- 
tively. In addition to fattening, so called Naik term and 
Lepage term are added. The Naik term is a product of 
three link variables along one direction 0. 

Uj? aik (x) = c N U^x)U^x + fiU^x + 2/x) (5) 

and the Lepage term is 

xUl{x + v + (i)Ul(x + jX). (6) 

We compare the action including staple term and 
the Naik term only (Staple+Naik action), i.e. the same 
as that of [8j, Asq action and Asqtad action. 
The Staple+Naik action is 



s a ^( x 'V) = E Vn(x)sign(ji) x [ 

x+fi 

+c N U^ k (x)6 y , x+3fl ] (7) 

where and c\ = 9/32, w 3 = 9/64 and cjv = —1/24. 

The Asq action contains S^ 3 \ S^ 5 \ and Lepage 
term and the Naik term but no tadpole renormalization 
factor. 

4 

s ap( x >y) = E Vi^{x)sign(fj,) x [ 

+ E w ? S $pt(x)) + c L U^(x))5 y , x+ , 

T 

+c N U^ k (x)S yM ] (8) 

where c x = 5/8, w 3 = 1/16, w 5 = 1/64, w 7 = 1/384, 
c L = -1/16 and c;v = -1/24. 

The full Asqtad action contains the tadpole renormal- 
ization factor. Different from we do not absorb one 
power of uq into the quark mass in our calculation of the 
quark propagator, and use c\ — 5/8itg , w 3 — l/Wu^ 3 , 
w 5 = l/64u^ 5 , w 7 = l/384u^ 7 , c L = -l/16u^ 5 and 
c N = -l/24u(7 3 . 



The inversion of the S~^(x, y) is performed via conju- 
gate gradient method after preconditioning 9] as shown 
in the appendix. 

We consider 

Pn{n,v) = \ E ^( a; ) 

H=l,--- ,4 

x {U* at {x)8 ytW+M + U» aik (x)5 yM 

- Ul at {x)^ y ^-U^ k {x)U y ^} (9) 

and solve for the bare valence quark mass mo, 

(p + m Q )<j ) = p (10) 

using the 3x3 color matrix expression to the source term 
P- 

The propagator can be expressed in the formpj 



s( q ) = zM q2 + M[q? 

Since tr^^ — 0, tr over color and flavor yields 
trS{ 9 ) = 16N, Z ^ )M[q) 



(11) 



q 2 + M(q) 2 
= 16N c B(q) 



(12) 



On the other hand 



t 



triiiS)=trZ 2[ q)- f+M{q)2 
q 2 + M{q) 2 



q* + M{qf 

7.Jn\ 

16AT C <7 2 
16N c q 2 A(q) 



The dynamical mass of the quark is 

B(q) 



M(q) 



A{q) 



and the quark wave function renormalization is 
A{q) 2 q 2 + B(q) 2 



Mi) 



A(q) 



(13) 



(14) 



(15) 



We remark that our method of deriving ^4(<7) and B(q) 
from the lattice data S(q) is different from [7J, however 
the results are equivalent. 



III. RENORMALIZATION EFFECTS 

In the continuum theory, the quark wavefunction 
renormalization Z^,{q 2 ) is defined by the renormalized 
inverse quark propagator S~ 1 (q) as [l(llllj 



S- l {q)=5 ab Z^{q 2 ){ii + M[q)) 



(16) 



On the lattice, renormalized operators must obey the 
same renormalization condition as the continuum theory. 
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The quark field renormalization factor is usually defined 
by the amputated Green function of the vector current. 

We define the colorless vector current vertex by 
using [Io| 



G„(q, P ) = / d^ye^y+^iqiyMx^qixMO)) 



and 



T fi (q,p) = S- 1 (q)G^q, P )S-\p + q) 



(17) 
(18) 



The general expression of the vertex of the vector current 
is 

r M (g) = <5 a , 6 {3i((7 2 )7 M + ig2(q 2 )q^ + gz(q 2 )q,j.i 

+m(q 2 )h„J}} (19) 

where for 16 flavors 
1 



(20) 



The Ward identity in the renormalized form tells 

. d 



(Tji)p(g) 



-S- R '(q) 



(21) 



After multiplying both sides by Z^, one obtains 

. d 



(22) 



where Zygi(q 2 ) = Z^(q 2 ) and Zy = 1 in the continuum 
limit and Z^{q 2 )/ gi {q 2 ) = Z^ IOM {a-\q 2 ) is expected 
to be independent of q in the limit of a~ 1 — > oo. 

In we observed that the fluctuation of the ghost 
propagator in the quenched and unquenched simulations 
differ significantly. 

The renormalization of the Green function for the 
ghost-gluon coupling is 



G{q,g) -> ~z^G(q,g'), 
?3{q,P,g) -» ziTa(q,p,g r ), 



fJ 



- -1/2 — 1 

ziz 3 z 3 g 



(23) 



When a fermion of mass m is coupled and the system 
is multiplicative renormalizable, the running coupling ex- 
pressed by the dimensionless structure function of gluon 
D and that of ghost D and the vertex function IV 

2 2 2 2 2 

-fP "I s / q P TYl 

- D ^^9) 1/2 xD(^,^g) (24) 

is invariant under the change of the cut-off A when g is 
transformed to 



In the Dyson-Schwinger equation(DSE) aproachp^ 
IbU l the running coupling g of the ghost-gluon coupling 



g(q) = Z^Zl /2 (v 2 , q 2 )Z,(n 2 , q 2 )g(n) (26) 

where Z\ is the ghost vertex renormalization factor, 
which is taken to be 1 in pQCD. On the lattice, how- 
ever, we took it as a scaling factor that adjust the r.h.s. 
of eq.® at the q = fx ~ 6GeV fits the value of pQCDQ. 
The value of 1/Zf is tabulated in TabkU 



TABLE I: The 1/Zf factor of the unquenched SU(3). 



config. 


heavy 


light 


N, 


comments 


MILC C 
MILC/ 


1.49(11) 
1.37(9) 


1.43(10) 
1.41(12) 


2+1 
2+1 


p mp = 6.83, 6.76 
P imp = 7.11,7.09 



We observed that running coupling g in the rhs needed 
to be enhanced by Z^ , so that it agrees with the pQCD. 

The running coupling of the quark gluon coupling is 
defined by the dimensionless structure function of gluon 
D and that of quark D^, and the vertex function T 3l p : 



2 2 2 2 2 

_p m ,q P m 

A 1 "' 5 ^ = 9 3V,( A2 ' A2' A2"' 5 ^ 

2 2 11 

n , q m n ,q m 



(27) 



When the quark propagator is not calculated consistently 
using the simplified action as compared to the real action 
used in producing the gauge fields, one should define the 
normalization of the quark propagator Z$ via other in- 
formation. 

Since at \x ~ 6GeV, the quark wave function renormal- 
ization factor in perturbative QCD (pQCD) is close to 1, 
we define the scale of Zjf,(q 2 ) as 



Z^q 2 ) 



Z vgi (p 2 )A(q) 
A(q) 2 q 2 + B(qy 



(28) 



where <?i(/i 2 ) is defined from the running coupling in the 
following argument. The function B{q) defines the quark 
inverse propagator Z^(q) and the dynamical mass M(q) 
as 



B(q) = A(q)M(q) cx 



M(q). 



(29) 



In the following, we define M(q) = A4(q) and renor- 
malize the quark field as ^/ZpiJ^j^bare — fpR- 

The Ward identity implies that the wavefunction 
renormalization can be defined via conserved vector cur- 
rent vertex. 

The running coupling of quark-gluon coupling is cal- 
culated as 



9 



1 1/2- 

z 1 z 3 ' z 3 g 



(25) 



7(9) 



Zf l zl /2 tf,<?)Z 2 ( M *,?)gOi). (30) 



4 



At 



the 



renormalization 



point q 



we 



fix 



Z 2 (fi 2 ,[i 2 ) — 1 and Z 3 (fi 2 ,ii 2 ) = 1, and thus 



Z X Z. 



-l 



(31) 



i.e. Z x = Z\ 

In [13, Zy is found to be almost constant for q > 
5GeV, and at fi ~ 6GeV we identify Z v gi(fi 2 ) with 
Zf(fj, 2 ). When quark wave function renormalization fac- 
tor is measured using the Asqtad action, the running cou- 
pling given by ea. ljSTH) is expected to agree with the pQCD 
result in high momentum region. But in the calculation 
with Staple+Naik action, the wavefunction renormaliza- 
tion is to be modified by the tadpole renormalization fac- 
tor uq. Further, in the approximate calculation, the Zf 
would be modified as Zi, and to give the consistent g(q), 
additional correction to Z 2 (q) would be necessary. 

We define Z 2 (q) = l/Z^q 2 ) = Z 2 {q) / 'Zy <?i(/j 2 ) and 
identify \jZyg\ (n 2 ) for the Staple+Naik action as the av- 
erage of the product of vertex renormalization factor and 
the tadpole renormalization factor i.e. Z x uq = 1/1.38 for 
MILC C and 1/1.36 for MILC/. 

Thus, the renormalization of Z$ on the lattice is de- 
fined by the renormalization of the running coupling on 
the lattice defined at fx ~ 6GeV. 

In the case of Asq action, the renormalization factor 
Z 2 (q) is suppressed by about 10% and the mass function 
M(q) is enhanced by about 10% as compared to the Sta- 
ple+Naik. By inclusion of the tadpole renormalization 
factor uq to the Asq action, the Z 2 (q) become consis- 
tent with those of Staple+Naik action with correction by 
Zygi([i 2 ). In 0, the scale of Z 2 (q) calculated by the 
Asqtad action is fixed to 1 at q = 3GeV, but we adopt 
the bare lattice data Z 2 {q) = Z 2 (q) and M(q) of Asqtad 
action to compare with those of Staple+Naik with correc- 
tion Zygi([i 2 ). Coincidence in the case of mo = 82.2MeV 
is shown in FigsQ]and|21 In the case of mo = 11.5MeV, 
the M(q) and Z 2 (q) of Staple+Naik in the infrared re- 
gion are shifted from those of Asqtad slightly, but the 
difference is about 10% at most, as shown in Figs. 3 and 
4, respectively. 
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FIG. 2: The Zi{q) of Staple+Naik action(stars) and the Asq- 
tad action(triangles) of MILC C with the bare quark mass 
m = 82. 2MeV. (Color online) 
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FIG. 3: The mass function M(q) of Staple+Naik action(stars) 
and the Asqtad action(triangles) of MILC C with the bare 
quark mass mo = 11. 5MeV. (Color online) 
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FIG. 1: The mass function M(q) of Staple+Naik action(stars) 
and the Asqtad action (triangles) of MILC C with the bare 
quark mass mo = 82. 2MeV. (Color online) 



FIG. 4: The Z 2 (q) of Staple+Naik action(stars) and the Asq- 
tad action(triangles) of MILC C with the bare quark mass 
m = 11.5MeV. (Color online) 
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IV. ANALYSIS OF M(q) AND Z^q 2 ) 



In perturbative QCD 
quark is expressed as@, 



CD), dynamical mass of a 



M(q) = 



47r 2 d M (ggUlog(<z7A 2 



QCD 



3g2[l g( M VA 2 CI5 )]^ 
m{^)[\og{^/A? QCD )]^ 



[\og{q>/kl CD )Y 



(32) 



where d,M 



I second term is the contri- 



33 - 27V/ 
bution of the massive quark. 

In this analysis of the lattice data, the quark con- 
densates — (qq(n)) and Aqcd are the fitting parameters. 
In the MILC/ lattice, the bare masses are 0.0062/a = 
13.6MeV and 0.0124/a = 27.2MeV for the u - d quarks 
and 0.031/a = 68.0MeV for the s-quark. In the MILC C 
lattice, the corresponfing masses are 0.007/a = 11.5McV 
and 0.040/a = 65.7MeV for the u — d quarks and 
0.050/a = 82.2MeV for the s-quark. 

The mass function ea. (|32|) is based on pQCD and can- 
not fit the data below 2GeV and we try the phenomeno- 
logical fit 0, [13 



M(q) = 



cA 3 



A 2 



mo 



(33) 



where Too is the bare quark mass. Parameters of c and 
A are summarized in TABLE [n] The fitting of the mass 
function of MILC C and that of MILC / are shown in Figs. 5 
and 6, respectively. 




FIG. 5: The dynamical mass of the MILC C quark 
(Staple+Naik) with bear mass mo — 11.5MeV(stars), 
65.7MeV(diamonds) and 82.2MeV(triangles) and the phe- 
nomenological fits. (Color online) 

We observe that as the bare quark mass becomes 
heavy, c becomes smaller but the product cA becomes 
larger. Although the MILC C configurations of bare mass 
too = 82.2MeV, with diffeernt (3 agree within errors, the 
MILC / configurations of bare mass too = 68MeV show 
dependence on (3. The mass function of j3 = 7.09 is 
smaller than that of (3 = 7.11. In the case of (3 — 7.11, 




FIG. 6: Same as FIG. but MILC/ quark mith bear 
mass mo = 13.6MeV(stars), 27.2MeV(diamonds) and 
68.0MeV(triangles) and the phenomenological fits. (Color on- 
line) 



TABLE II: The parameters c and A. 





mo(MeV) 


c 


A(GeV) 


cA(GeV) 


6.76 


11.5 


0.44(1) 


0.87(2) 


0.383 




82.2 


0.30(1) 


1.45(2) 


0.431 


6.83 


65.7 


0.33(1) 


1.28(2) 


0.420 




82.2 


0.30(1) 


1.45(2) 


0.431 


7.09 


13.6 


0.45(1) 


0.82(2) 


0.368 




68.0 


0.30(1) 


1.27(4) 


0.381 


7.11 


27.2 


0.43(1) 


0.89(2) 


0.383 




68.0 


0.32(1) 


1.23(2) 


0.397 



the chiral limit M(0) is consistent with that of MILC C 
and we find M(0)=0.37(l)GeV. However, in the case of 
/3 = 7.09, too = 13.6MeV, the lowest three momentum 
points of M (q) are systematically smaller than the other 
points. The slope of the (3 = 7.09 and that of 7.11 are 
almost the same. A preliminary analysis of cA using the 
Asqtad action suggests that the data of (3 = 7.09 and 
MILC C obtained by the Staple+Naik are underestima- 
tion, and the discrepancy between MILC C and MILC/ 
of about 0.02GeV remain. We expect that the M(0) in 
the continuum limit would be about 0.38GeV, consistent 
with the value obtained by meromorphic parametriza- 
tion of lattice data^^. In an DSE approach p^j. devia- 
tion of the mass function M(0) from a linear function of 
to(C = 19GeV) where £ defines the scale of the system 
is claimed. Although to(£ = 19GeV) is not identical to 
to , M(0) — m(( = 19GeV) of the DSE increases near 
the chiral limit in contrast to our naive fitting shown 
in FIG[7| which implies that the chiral symmetry effect 
monotonically decreases as the bare mass increases. 

We show the lattice results of ^2(9) of MILC/ and 
MILC C in Figs. |Hland|51 respectively. 

The apparent difference in the formulae of Q and our 
work are only in the expression and in fact they are equiv- 
alent. The Z.2{q) agree with each other. 
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FIG. 7: The chiral symmetry breaking mass cA = M(0) — mo 
as a function of bare mass and its chiral limit. Dotted line is 
the extrapolation of MILC/ (3 — 7.09, dashed line is /3 = 7.11 
and the dash-dotted line is that of MILC C . (Color online) 
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FIG. 8: The Z 2 (q) of MILC/ with bear mass m = 
13.6MeV(stars), 27.2MeV(diamonds) and 68MeV (triangles). 
(Color online) 



The A(q) denned as S(q) = 



1 



A{q)i-B{q) 

parametrized in the operator product expansion ap 



is 
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FIG. 9: Same as FIG. H but quark of MILC C with 
bear mass mo = 11.5MeV(stars), 65.7MeV(diamonds) and 
82.2MeV(triangles). The last two almost overlap. (Color on- 
line) 



proach as 2 

Mi) 



= i 



N c q 2 



7ra s (/u 2 )(F 2 ) M 37ra s (/x 2 )(g 5 4(7) A 



3N c q 4 



(34) 



where (F 2 )^ is the gluon condensate, a s (fj, 2 )(A 2 ) ll is the 
A 2 condensate and (qgAq)^ is the mixed condensate, 
which are the fitting parameters. 

Similar parametrization was done by 2] , in which /i = 
2GeV and A(q) for q > lGeV is parametrized as 



(35) 



In this parametrization, we obtained ci = 0.25GcV 2 , 
c 2 = 0.061GeV 4 in the case of mn = 27MeV. Comparing 
with data of quenched simulation 2\{c\ = 0.37±0.06GeV 2 
and C2 = — 0.25 ± 0.04GeV 4 ), ci is the same order but c 2 
is smaller and have opposite sign. The result depends on 
the choice of the renormalization point [i. 

The pQCD result of the wave function renormalization 
factor Z 2 (q) in MS scheme up to four loop is given by 
12 2L |23|. Orsay group expanded the af s in terms of 



„MOM 



(?) 



(vi- 



and obtained Z^ ert (q" 



1/Z 2 (q) as a 



function of a and fitted the Z^(q 2 ) of Wilson overlap 
fermion in the MOM scheme 

{Z^\u?))^S-\q) = (^V))- 1 ^?) 



d(q 2 /n 2 ,a(n)) 



A{N 2 



ii>5 a b 



' ^ £r V) , d{g 2 /n\a{»)) (A 2 ), 

rypert 



(m 2 ) 



4(iV c 2 - 1) 



+ ■ 



(36) 



where d(q 2 / ji 2 , a(/i)) is the solution of the renormaliza- 
tion group equation 



{(-70+7V) 



47T (ifog/i 2 

which can be written as 



}d(q 2 /fi 2 ,a(fi)) = (37) 



d(q 2 /v 



2 a( M )) = d(l,a(?))' " (/ " 



(-T0+742)/^0 



The anomalous dimension of A 2 in the lowest order 

is[H0[23 7 i o 2 = !-! 

(A 2 ) as a fitting parameter and calculate 



a{q) 

f A 2 
nf and 7q 



(38) 



0. We adopt 



Z^(q 2 ) 



Zvgiiv 2 



qa\ (-70+7 A 2)//3o 



ZT\q 2 ) 



(A 2 ), 



A(N 2 - 1) 



7 pert 



£2 

rA 



(M 2 ) 
(39) 
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where a(q) are data calculated in the MOM scheme us- 
ing the same MILC/ gauge configuration^. 

We fitted the Z^q 2 ) of MILC/, to = 27.2MeV data 
and choosing /i = 2GeV and = 0, obtained (A 2 )^ = 
1.6(3)GeV 2 , which is compatible with the Orsay group 
data 2.4 ± 0.3GeV 2 for the Wilson fermion. 

In FIG I1UI we show the lattice data of 

Z^,(q 2 ) = ttt~T ~ from the Staple+Naik action 

with Zvgi{p 2 )=^-2>& i.e. the average of the running cou- 
pling renormalization factor of MILC/ in the previous 
analysis^, and the phenomenological fit Z^ lt (q 2 ) by 

eq.dSHl- The pQCD result Z v * Ti \q 2 ) above q = lGeV 
is shown by the dashed line and the lattice data of 
a s (q) below q — 0.5GeV obtained by the ghost-gluon 
coupling arc also plotted. The difference between a s (q) 
and Z^{q 2 ) in the infrared implies that Zy is strongly 
momentum dependent in the infrared [icij. 
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FIG. 10: The quark wave function renormalization Z^(q 2 ) 
of MILC/ bear mass mo = 27.2MeV. The solid line is the 
phenomenological fit as [Tfll | and the dashed line is the pQCD 
result. The points below 0.6GeV are gi(q 2 ) approximated 
by the running coupling obtained by the ghost-gluon vertex. 
(Color online) 

The suppression of Z^(q 2 ) in the infrared suggested by 
the behavior oiaJq) could be an artefact. In the case of 
Wilson fermion 26], the Z^(q 2 ) in the continuum limit is 
infrared finite. 



V. CONCLUSION AND DISCUSSION 

We measured the quark propagator of the MILC / lat- 
tice (28 3 x 96) of /3 lmp = 7.09 and 7.11 and MILC C lattice 
(20 3 x 64) of p mp = 6.76 and 6.83 using the Staple+Naik 
action. After the renormalization using the information 
of the running coupling and the tadpole renormalization, 
we obtained a good agreement with the lattice bare data 
of the Asqtad action in the case of MILC C . 

We observed that the denominator of the quark renor- 
malization factor 



A(q) 2 q 2 +B(qY 
A(q) 



is a steeply increasing function as q approaches which 
causes the infrared suppression. Although the function 
B(q) in the numerator is also an increasing function of 
q as q approaches 0, it is almost constant as the bare 
mass ttiq decreases, in contrast to A(q) which becomes 
larger as the bare mass too decreases. In the DS equation 
approach [Tsj. a v (q) and <J s (q) correspond to our A(q) and 



B(q), respectively. cr v (q) corresponds to 



and 



q* + tyIq 

its mo dependence is qualitatively consistent with our 
lattice data. The dominant term of cr s (g) corresponds 
cA 3 

to — — n . , „ r"oT an d its too dependence is weak in 

(q 2 +m 2 )(q 2 + A 2 ) 

the q ~ O.lGeV region consistent with our lattice data. 
The infrared feature of Z^,(q 2 ) = — , ; is related to 



the Kugo-Ojima color confinement criterion: 



1 



Zz 



Zi 
Zz 



Zi-tp 

~z^ 



= 



(41) 



(40) 



i.e. the Kugo-Ojima parameter u = — 1. 

In the quark sector, divergence of Zi(tq) or vanishing 
of Z 1 p(q 2 ) in the infrared is consistent with the Kugo- 
Ojima confinement criterion. Infrared finiteness of Z2 
implies that is infrared vanishing, and infrared finite 
Z^p supports the idea that the running coupling freezes 
to a finite value in the infrared. This situation is similar 
to that of the gluon sector. The infrared finite lattice 
data of Z3 implies that Z\ is also infrared vanishing. In 
this case running coupling defined by H26JI vanishes un- 
less Z3 is infrared vanishing. Definition of the running 
coupling from the triple gluon vertex is more ambiguous 
than the definition (|26[) from the ghost-gluon vertex, and 
we expect that the fluctuation of the ghost propagator in 
the infrared induces the artefacts in the infrared running 
coupling. 

We compared the Staple+Naik action and the Asqtad 
action and found that the mass function and the quark 
propagator are consistent with each other. They are con- 
sistent with the results of other group[3 0,113. 

We confirmed A 2 condensates which was found by the 
Orsay group in the running coupling of quenched lat- 
tice simulation [281 an d observed in the unquenched lat- 
tice simulation[l|, [l(| also in the quark propagator. We 
observed also the indication of the qq condensate. 

The effect of dynamical chiral symmetry breaking is 
seen in the mass function M(q). By extrapolation to 
the chiral limit we obtain M(0) ~ 0.38GeV, which is 
about 20% larger than Q, but consistent with DS equa- 
tion approach|l8j. An origin of the difference from 
is their renormalization of the mass function M(q) at 
3GeV to a theoretical value, which we did not adopt. In 
the definition of the quark mass from the lattice data, 
one should take into account the mass renormalization 
by taking into account the lattice constant a dependence 
of the propagator |H HI HJ EJ HI • 

Calculation of the quark propagator of MILC / using 
the Asqtad action is underway and the detailed compar- 
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ison of the Asqtad action data and the DS equation ap- 
proach will be published in the forthcoming paper. 
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APPENDIX A: THE LATTICE CALCULATION 
OF THE KS FERMION PROPAGATOR 

The Dirac gamma function of KS fermions are defined 
asQ 



where 



0' 



(/*) _ J 1 v < /* 



otherwise 



(1) 



(2) 



In momentum space the inverse propagator is expressed 
as 



S a p(Pi m ) = i ^2(lf^)ai3{^smp^ - ^-sinSp^) +m8 af} 

(3) 



where 



tap = n s ait0ll 



(4) 



The algorithm for calculating the propagator is as fol- 
lows. We define the operator devided by the bare mass 
of the dynamical quark m as 



M=(I+-p) 
m 

and decompose even sites and odd sites, 



(5) 



M 



I ij 



I -L-U (6) 



where 



L = 



and 



U = 





-±Peo 



-hPoe 



(7) 



v 

Using the Eisenstat trick, we define 



(8) 



M = {I - L)- l M(I -U)- 1 = (I + L)(I-U-L)(I + U) 
', \ 

i-(±fp eo p oe 



I-LU 



(9) 



We note that M is hermitian and the conjugate gradi- 
ent method and/or BiCGstab method are applicable. 
With use of the definition 



f 

m 



(10) 



and 



we solve 



(I 



-Wo 



(11) 



(12) 



= P, 



The odd site solution is 



with initial value 

4>o = P' - _ 

The decomposition into even sites and odd sites is done 
by multiplying the projection operator P e and P a , and 
the color source p' is taken as a vector of 3 x 3 matrices. 

The convergence condition for the conjugate gradient 
iterations is 



\\p\\ 



< a few per cent at most 



(13) 



where the used norm is maximum norm in the space of 
site, color and flavor, and the accuracy gets 10 _1 higher 
if L 2 norm is used. 
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